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ABSTRACT 

We model the mass distribution in the recently discovered Einstein ring LBG J213512.73- 
010143 (the 'Cosmic Eye') using archival Hubble Space Telescope imaging. We reconstruct the 
mass density profile of the z — 0.73 lens and the surface brightness distribution of the z = 3.07 
source and find that the observed ring is best fit with a dual-component lens model consisting 
of a baryonic Sersic component nested within a dark matter halo. The dark matter halo has an 
inner slope of 1.42±g;||, consistent with CDM simulations after allowing for baryon contraction. 
The baryonic component has a mass-to-light ratio of 1.71 jlQ'|| Mq/Lbq which when evolved to 
the present day is in agreement with local ellipticals. Within the Einstein radius of 0.77" (5.6 
kpc), the baryons account for (46 ± 11)% of the total lens mass. External shear from a nearby 
foreground cluster is accurately predicted by the model. The reconstructed surface brightness 
distribution in the source plane clearly shows two peaks. Through a generalisation of our lens 
inversion method, we conclude that the redshifts of both peaks are consistent with each other, 
suggesting that we are seeing structure within a single galaxy. 



Subject headings: cosmology: observations — cosmology: dark matter 
ular — galaxies: individual: LBG J213512. 73-010143 



galaxies: elliptical and lentic- 



1. Introduction 

The measurement of galaxy mass profiles has 
proven a powerful observational probe for testing 
models of structure formation. In recent years, 
the slope of the inner mass profile has become a 
subject of much contention with the results of cold 
dark matter (CDM) simulations being discrepant 
with observations. At the present time, this dis- 
crepancy persists. 

Navarro. Frenk fc White I (|l996h first proposed 



an analytic approximation to describe the mass 
density profiles of halos in their CDM simulations. 
At a radius much smaller than a particular scale 
radius, the mass density of this so called 'NFW 
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pr ofile scales as p(r) oc r . Later simulations 
by iMoore et al. ( 1998 . 19991 ) indicated a steeper 
inner slope. This gave rise to the 'generalised 
NFW (gNFW) profile which at radii smaller than 
the scale radius follows p(r) oc r~ a with values 
of a around 1.4 to 1.5. The most recent simu- 
lations boasting a significantly higher resolution 
now conver ge on a slope somewhere in the range 
UL>a<1.2 ([Navarro et "ai~ll2004l:lDiemand et al 



2005) 



These results appear to strongly disagree with 
measurements of the inner slope from dynamical 
studies. Several groups measuring rotation curves 
of low surface brightness galaxies (LSBs; these are 
believed to have a high dark matter fraction and 
therefore be minimally affec ted by baryons) fi nd 



a range of slopes; < a < 1 (|de Blok et al. 2001 
Ide Blok fc Bosnia I l200a ISwaters et al. I boo 



)3; 
1. 



Sock kens fc Giovanelli I |2005J). lHavashi et al 
(|2004h purported that this discrepancy could be 
reconciled by directly comparing against rotation 
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curves of simulated halos. Howev er, this was 
strongly rejected by Ide Blok I (j2005h who found 
that only one quarter o f the 5 1 galaxies in the 



study of lHavashi et al. ( 2004 ) were consistent 



with CDM . In the latest episode of this ongo- 
ing debate, Havashi fc Navarro ( 20061 ) claim that 
non-circular motions in simulated CDM halos aris- 
ing from a triaxial potential can explain the range 
of measured rotation curves seen in LSBs. Clearly 
dynamical measurements, whilst potentially very 
powerful, have many complexities preventing a 
straight-forward interpretation. 

Gravitational lensing has for some time now 
provided an attractive alternative means of mea- 
suring mass profiles without the difficulties asso- 
ciated with dynamical techniques. Primarily, this 
is motivated by the simple fact that the deflection 
angle of a photon passing a massive object is inde- 
pendent of the dynamical state of the mass within 
the object. Strong lensing systems with multi- 
ple images of a background source can constrain 
the radial profile of the projected mass density 
of the lens by searching for the best fit to the ob- 
ser ved image positions (for example, see the review 
by ISchneider. Kochanek fe Wambsganss" 20061). 



This t echnique was enhanced by I Sand. Treu fe Ellis 



(|2002l ) to incorporate extra constraints from the 



observed velocity dispersion profile of the lens and 



dTreu & KooDman 


s 2002: Koopmans & Treu 


20031 Sand et al. 


20043 . 


Dve & Warren 


(2005, hereafter DW05) showed 


how Einstein ring systems, i.e., strong lens sys- 



terns where an extended source is imaged into a 
complete or near-complete ring, can place stronger 
constraints on the mass profile of the lens than sys- 
tems with multiple point-like im ages. This work 



used the semi-linear method of IWarren fc Dye 
(|2003l ) which has als o been used by several 



other studies to date dTreu fc Koopmans 



Treu et al. I 120061: iKoopmans et al. I 120061 ). A 



2004; 



Bayesian versi on of the semi-line ar method was 
developed by Suvu et al. (12006 ) and later en- 



hanced bv lBarnabe fc Koopmans 



( 2007t ) with the 



inclusion of linear constraints from stellar velocity 
moments. 

In this paper, we follow the procedure of DW05. 
We apply the semi-linear method to reconstruct 
the lens mass profile and source surface bright- 
ness image of the extraordinary Einstein ring 



system LBG J2135 -0102 recently discovered by 
Smail et al. ( 2007L see also Coppin et al. 2007). 
Due to the lens amplification, this system is of par- 
ticular interest as it represents one of the brightest 
examples of Lyman break galaxies known. 

We compare the fit given by five different lens 
models. With our most general model which com- 
prises a gNFW halo that hosts a mass component 
following the lens galaxy light, we separate the 
baryonic and dark matter contributions to the pro- 
jected lens mass within the Einstein radius. We 
also describe a simple extension of the semi-linear 
method enabling reconstruction of sources in mul- 
tiple planes at different redshifts. Our analysis 
adds to the growing list of Einstein ring systems 
that are now beginning to provide strong con- 
straints on CDM simulations. 

The layout of this paper is as follows. In the 
following section we briefly describe the data. Sec- 
tion[3]outlines the semi-linear method and our lens 
models. We present the results of our source and 
lens reconstruction in Section 01 Finally, we con- 
clude with a summary and discussion in Section[5] 
Throughout this paper, we assume the following 
cosmological parameters; Ho = 70kms _1 Mpc -1 , 
tt m = 0.3, fl A = 0.7. 

2. Data 

The image data were acquired as part of the 
Snapshot programme GO #10491 (PI: H. Ebel- 
ing) on May 8th 2006 using the Hubble Space 
Telescope (HST) Advanced Camera for Surveys 
(ACS). Imaging of the Einstein ring itself was 
completely serendipitous, the intended target be- 
ing the cluster MACS J2135. 2-0102 lying approx- 
imately 75" away to the south. The top panel of 
Figure [T] shows a section of the reduced ACS im- 
age with the Einstein ring and lens galaxy just to 
the left of centre. The cluster centre lies at 24° in 
a counter-clockwise direction measured from the 
+ve cc-axis. The nearby galaxies labeled A, B and 
C are likely cluster members. We consider the per- 
turbative effect from these nearby galaxies and the 
large scale shear from the cluster itself on the lens 
solution in Section 14.11 

The image is a stack of three dithered 400s ex- 
posures taken with the ACS Wide Field Camera 
(WFC) in the F606W filter. The stack was pro- 
duced using the STScI MULTIDRIZZLE pack- 
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age (V.2.7) thereby correcting for the geometric 
shear introduced by the WFC and giving a result - 
ing pixel sca l e of .05". The reader is referred to 
Smail et al. (2007) for more details. 
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Fig. 1— The ACS Wide Field Camera data ob- 
served in F606W. North lies at an angle of 197.4° 
measured counter-clockwise from the positive x- 
axis. Top: The Einstein ring and lens galaxy show- 
ing the nearby likely cluster members labeled A, 
B and C. The cluster MACS J2135. 2-0102 lies 75" 
away to the south as indicated. Bottom: The ring 
after subtraction of the lens galaxy satisfactorily 
fit with a Sersic profile. 



Keck spectroscopy taken on June 30th 2006 and 
July 24th 2006 unambiguously identifies the lensed 
source as a Lyman break galaxy lying at a redshift 
of z = 3.07. Similarly, the lens galaxy spectrum 
exhibits several strong absorption and emission 
features placing it at a redshift of z = 0.73. The 
total magnitude of the lens galaxy i n the F606W 
filter is Rfgogw = 22.34 ± 0.15 (|Smail et al 



2007 ). In additio n to the ACS imaging data, 



Smail et al. I (|20070 obtained a K' image using the 
Near Infra-Red Camera at Keck II on July 4th 
2006. The {Rfgogw - K')=2.8 colour of the lens 
and the spectral continuum is consistent with that 
of an early type spiral or SO galaxy. 

Before reconstructing the source surface bright- 
ness distribution, any contribution of flux from the 
lens galaxy to the ring flux must be removed. Al- 
though the lens contribution in this case is unlikely 
to be of concern, we carried out this exercise any- 
way to establish the visible morphology of the lens 
for our modelling later. An added complication in 
this case is that the lens galaxy light is sheared by 
the foreground cluster and this must be accounted 
for in the lens model (Section 13.3.11) . We fitted an 
elliptical Sersic profile ( Sersic T 1968) of the form 



L = L 1/2 cxp{-B(n)[(r/r s ) 1 /"-l]} (1) 

to the lens galaxy having first masked out pix- 
els within an elliptical annulus containing the ring 
(see Figure[3]). For the fitting of this profile and for 
the modelling later, we created a noise map based 
on the CCD gain, read noise and photon count. 
The parameters "L\/2, n and r s were allowed to 
vary in the fit as well as the elongation (i.e., major 
axis divided by minor axis), e, and the orientation 
(angle in degrees of the major axis to the pos- 
itive image cc-axis in a counter-clockwise sense), 
9. In the fitting, we convolved each trial sur- 
face brightness profile with the WFC point spread 
function (PSF) mode lled by the TinyTim package 
(|Krist fc Hook ll2004f ) for the F606W filteiQ. The 
best fit was achieved with the parameters listed in 
Table [1] which give an acceptable x 2 . 

We work in terms of the rest-frame B band lu- 
minosity in this paper for ease of comparison with 
other studies. The conversion from F606W flux 
(approximately rest-frame U) to B band luminos- 



1 The TinyTim PSF includes CCD pixel charge smearing and 
geometric distortion of the ACS image plane. 
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ity density was made using a k-correction and a 
colour correction (see Section l373|) . The fitted pro- 
file was subtracted from the unmasked image and 
the result is shown in the lower panel of Figure [1] 

3. Method of Analysis 
3.1. Semi-linear inversion 

For a full description of the sem i-linear inver- 
sion m ethod, the reader is referred to lWarren fc Dye 
(|2003l ) and DW05. We provide a very brief outline 
here for completeness. 

The inversion assumes the source plane and the 
image plane are pixelised. The manner in which 
the source plane is pixelised is not restricted, al- 
lowing the 'adaptive gridding' described in DW05 
whereby smaller pixels are concentrated in regions 
where there are stronger constraints (see next sec- 
tion). 

For a given lens model parameterisation, a PSF- 
smeared image is computed for every source pixel. 
All images are created using unit surface bright- 
ness source pixels. The problem of finding the 
factor required to scale each image such that their 
co-addition best fits the observed image is a linear 
one. These scale factors are the best fit surface 
brightnesses of the source pixels for the given lens 
model. The solution is written as 

s = (F + AH)~ 1 c, (2) 

where the square matrix F and the vector c have 



the elements 

Fife = ^ fijfkj/v 



J2fv d iH ( 3 ) 



Parameter 


Minimised Value 


Ll/2 


(3.14 ±0.25) x 10 U L SO /D" 


n 


2.90 ±0.12 


r s 


0.57" ± 0.09" 





110.1° ± 2.6° 


e 


1.38 ± 0.07 



Table 1: Parameters of the Sersic profile fit to 
the lens galaxy. Reading from top to bottom, 
these are; the surface brightness at the scale ra- 
dius (Li/ 2 ), the Sersic exponent (n), the scale ra- 
dius (r s ), the orientation of the major axis to the 
positive image x-axis in a counter-clockwise sense 
(9) and the elongation, i.e. major axis divided by 
minor axis (e). 



and s is a vector containing the best fit source pixel 
surface brightnesses. Here, dj is the observed flux 
in image pixel j, Oj its error and f%j is the flux 
in pixel j of the image of source pixel i for the 
current lens model. The solution is regularised 
by the square regularisation m atrix H, scaled by 
the regularisation weig ht A fsee lPress et al. II200U 
and Warren & Dye 2003). The standard errors of 
the reconstructed source pixels are given by the 
diagonal terms of the covariance matrix which is 
just: 

(F±AH)- 1 . (4) 

This linear inversion is nested inside an outer 
loop that varies the non-linear lens model param- 
eters. For each trial set of lens parameters, the 
image of the reconstructed source is compared to 
the observed image and the x 2 computed within 
an elliptical mask containing the ring. The mask 
is designed to ensure that it includes the image of 
the entire source plane, with minimal extraneous 
sky. This means that only significant image pixels 
are used in the fit, making x 2 more sensitive to 
the model parameters. 

The regularisation scheme is chosen to best suit 
the reconstructed sour ce surface brightness image. 
Warren fc Dye I (|2003l ) showed using simulations 
that linear regularisation works well for this appli- 
cation. However, as we discuss in the next section, 
regularisation biases the solution between different 
lens models and so we set A = and use an adap- 
tive source grid when searching for the best fit lens 
parameters. 

In DW05, we used a downhill simplex method 
to find the set of lens parameters that minimises 
X 2 . However, the v 2 surface is pitted with lo- 



cal minima and noisy. IWarren fc Dvel (|2003ft de- 
scribed how this noise can be reduced by sub- 
gridding the image pixels when computing the 
source pixel images. In this work, we therefore 
divide image pixels into 16 x 16 sub-pixels. To 
prevent becoming trapped in local minima, we use 
a simulated annealing downhill simplex minimisa- 
tion algorithm. We find that a slow exponentially 
cooled temperature with a half-life of ~ 50 iter- 
ations works extremely well in finding the global 
minimum. 
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3.1.1. Multiple source planes 



3.2. Adaptive source plane grid 



As we show in Section I4.2[ the reconstructed 
source surface brightness distribution of LBG 
J2135-0102 has two distinct peaks. An inter- 
esting question to ask is therefore whether these 
peaks are two separate objects at different red- 
shifts or, given their small angular separation, 
two flux peaks withi n the same object. The ring 
spectra presented in Smail et al. I ( 2007 ) show no 
evidence of a different secondary source redshift 
although the fainter source would be dominated 
by the brighter and much more highly magnified 
primary source (see also the discussion of the CO 
emission from this system in lCoppin et al. Il2007[ ). 
We explain here how the semi-linear inversion 
method can be modified in a simple way to place 
constraints on this problem. 

The inversion described above is completely 
general in terms of the geometry of each source 
pixel and the lens model. Qualitatively, this means 
that each source pixel image can be created as- 
suming a different lens configuration. The practi- 
cal use of this is not in changing the lens mass 
profile, but the lens geometry so that different 
source pixels can be assigned different redshifts. 
In the case of two source planes, this quantita- 
tively means that the vector s in equation (|2|) 
holds two groups of source pixels, s\ — > s n and 
Sn+i — * Sn+m- The matrices F and H and the 
vector c are c omputed in exactly th e same way 
as described in Warren fc Dye ( 20031 ). except the 
lens configuration for the first n pixels uses a dif- 
ferent source redshift to that for the remaining m 
pixels. 

In this paper, we carry out the usual recon- 
struction with a single source plane, then a sepa- 
rate analysis with two source planes. For the sec- 
ond analysis, we define two source planes; plane 
1 containing pixels that belong to the brightest 
and more highly magnified source and plane 2 
containing pixels of the fainter source as Figure 
[3] shows. To investigate whether the redshift of 
source 2 is consistent with that of source 1, we 
hold plane 1 fixed at the measured spectroscopic 
redshift z = 3.07 and allow the redshift of plane 2 
to vary. This introduces an extra non-linear lens 
model parameter (see Section I3.3[) . In the min- 
imisation, all lens parameters including this addi- 
tional parameter are allowed to vary. 



We adopt the gridding algorithm of DW05 in 
the current paper which creates smaller source 
plane pixels where the magnification is higher. 
The gridding algorithm starts with a regular mesh 
of large pixels. For a given lens model, the av- 
erage magnification /i^ of every source pixel i is 
computed. Those pixels that meet the criterion 
fJ'ifi > s are then split into quarters, where r$ is 
the ratio of the area of pixel i to the area of an 
image pixel and s is the 'splitting factor'. Hav- 
ing finished the initial loop through all pixels, the 
process is repeated for the sub-pixels, then for the 
sub-sub-pixels and so on until all pixels satisfy the 
splitting criterion. The entire procedure (which 
takes a fraction of a second on a modern desktop 
computer) is repeated every time the lens model 
changes. 

The splitting factor is set empirically by mea- 
suring the improvement in the fit brought about 
by the splitting for different values of s with a 
fixed lens model close to the best solution. The 
splitting factor is successively reduced until no sig- 
nificant improvement is obtained. For our unreg- 
ularised reconstruction, we find a splitting factor 
of s = 9 satisfies this condition. By initialising 
the adaptive gridding with a source plane of regu- 
lar 0.1" x 0.1" pixels, this splitting factor ensures 
that images of split source pixels do not exceed the 
Nyquist sampling resolution set by the image PSF 
(see DW05 for more details on how s is set and how 
the initial grid size is chosen). This is important 
since it keeps degeneracies between source pixels 
to a very low level, thereby allowing the number 
of degrees of freedom to be accurately determined 
(see below). 

The main advantage of an adaptive grid is that 
the resolution of the reconstructed source can be 
enhanced relative to a regular grid without the 
need for regularisation. Regularisation has the 
effect of smoothing the reconstructed source light 
profile, reducing the effective total number of pa- 
rameters and hence increasing the number of de- 
grees of freedom, by an amount that cannot be 
satisfactorily quantified. This is especially prob- 
lematic when comparing different lens models, as 
a fixed regularisation weight for one model gener- 
ally does not give the same increase in number of 
degrees of freedom for another. This problem was 
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noted by iKochanek, Schneider fc Wambsganss 
(|2004h . Adaptively sized pixels extract maximum 
information from the lens image without need of 
regularisation. Therefore, provided the degenera- 
cies between all pairs of source pixels are negligi- 
ble, the number of degrees of freedom of the fit 
is a well-defined number. This allows, firstly, un- 
ambiguous assessment of whether a given model 
provides a satisfactory fit to the data and secondly, 
unbiased comparison of different model fits. 

Suyu et al. ( 20061 ) and Barnabe fc Koopmans 



(|20071 ) suggest an alternative means of circumvent- 
ing this problem with a Bayesian approach. They 
point out that the Bayesian evidence enables dif- 
ferent lens models and regularisation weights to be 
compared in an unbiased way. However, we prefer 
to adhere to the adaptive grid technique without 
regularisation since this has the appealing charac- 
teristic that the error image of the reconstructed 
source is much more uniform. To ensure that the 
number of degrees of freedom we compute in our 
unrcgularised solution is accurate, we verify that 
the degeneracies between all source pixel pairs are 
negligible by ensuring that all off diagonal terms 
in the covariance matrix given by equation Q are 
negligible. 

Although regularisation is not used when quan- 
tifying lens model fits, we use it in this paper solely 
for the aesthetic purpose of obtaining a higher 
resolution source. In DW05, we used zeroth or- 
der regularisation to simplify the implementation. 
However, this is an over-simplistic type of regular- 
isation since it poorly represents real astronom- 
ical sources. A better scheme is linear regular- 
isation where pixels that differ more strongly in 
surface brightness from their neighbours are pe- 
nalised more heavily. We implement a form of 
linear regularisation on the adaptive grid in the 
current paper by constructing a regularisation ma- 
trix that takes the difference between a given pixel 
i and the sum of all neighbouring pixels j weighted 
by Wij = (aj/a,i) Aexp(— y^/2a 2 ). Here, a is 
the source pixel area, yij is the separation of the 
centres of pixels i and j and N is a normalisa- 
tion constant set such that 



1. We 

■ — •JU-T" 

set a — 0.05 . As described in DW05, we use 
a smaller splitting factor when applying regular- 
isation. Figure [3] shows the adaptive grids upon 
which the regularised and unregularised source is 
reconstructed. 



3.3. Lens model 

We compare the fit given by five different pa- 
rameterisations of the lens mass profile. The most 
general of these is a dual component model. In this 
model, one component follows the lens galaxy light 
(the baryons) and the other, a gNFW profile rep- 
resents the dark matter halo. The remaining four 
single component models are special cases of the 
dual component model; 1) a singular isothermal 
ellipsoid (SIE), 2) a gNFW profile, 3) a power-law 
model, 4) a mass-to- light (M/L) model where the 
mass follows the lens galaxy light. 

The dual component model is that used by 
DW05. For the baryonic component, only the B 
band M/L ratio, is allowed to vary during the 
minimisation. Its morphology is fixed by the lens 
galaxy surface brightness profile that would have 
been observed in the absence of cluster shear. This 
is simply a transformation of the Sersic profile fit 
to the observed surface brightness profile (see Sec- 
tion [3XT]). 

The dark matter halo component is described 
by an elliptical gNFW profile with a volume mass 
density given by 



p(r) =Ps(l + r/r ) a - 3 (r/roy 



(5) 



The halo has seven parameters; the normalisa- 
tion, p s , the inner slope a, the offset of the halo 
centre from the baryonic centre, (Axh, At/ft), the 
scale radius ro, the elongation eh (equal to the 
major:minor axis ratio) and the orientation of 
the major axis to the positive image x-axis in a 
counter-clockwise sense, Oh- As we discussed in 
DW05, the scale radius is very weakly constrained 
in this model, hence we hold it fixed at the value 
of r n = 3.0 " in accordance with the results of 



Bullock et al. I (|200ll ) for a galaxy with similar to- 
tal mass and redshift properties to the lens. 

We include the effect of the nearby cluster in 
two ways. Firstly, we incorporate an external 
shear that is allowed to vary in magnitude, 7, and 
direction, d 1 (measured counter-clockwise from 
the +ve x-axis). The cluster convergence must 
also be accounted for, as this causes magnification 
of the ring image. We therefore calculate a con- 
vergence based on 7 and # 7 assuming a SIE pro- 
file (using a NFW profile for this purpose instead 
makes little difference to the best fit lens model 
parameters compared to the size of their errors). 
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Secondly, we model perturbations in the cluster 
potential from the nearby cluster galaxies. In this 
case, each galaxy is modelled as a SIE with a mass 
given by its total integrated light multiplied by a 
common B band M/L, \I/ CS . 

In summary, our dual component model has a 
total often free parameters; "J, p s , a, (Axh, Ay^), 
Oh, Gh, 7 , 0-y and ^ cg . The single component mod- 
els are obtained by holding various combinations 
of these parameters fixed. For the gNFW model, 
we only fix ^ = and allow all other parameters, 
including the cluster parameters (7, # 7 and ^ cg ) to 
vary. Similarly, for the power-law and SIE models, 
the halo scale radius, ro, is held at an arbitrarily 
large value having fixed = with the additional 
constraint a = 2 for the SIE. The pure M/L model 
allows only \& and the cluster parameters to vary. 

The deflection angles of both the Sersic and 
gNFW profiles must be numerically evaluated. 
Also, the surface mass density of the gNFW pro- 
file must be obtained by numerically integrating 
along the line of sight. We showed in DW05 how 
we convert the circular Sersic and gNFW profiles 
to elliptical profiles. 

The surface mass density of the Sersic profile 
is obtained by multiplying the fitted light pro- 
file given in equation ([I]) by the M/L ratio, 
To convert the lens galaxy flux measured in the 
F606W filter to a B band luminosity, we com- 
puted a k-correction and a colour correction us- 
ing a standard SO and Sa ga laxy template taken 
from iMannucci et al. I (|200lh - These SEDs are 
consistent wi t h the (R — K') colour measured by 



Smail et al. 



( 200 it and the observed lens mor- 
phology. The k-correction and colour correction 
for the SO template are respectively -1.85 mag and 
2.03 mag, and for the Sa template, -1.69 mag and 
1.97 mag. We take the average k-correction and 
colour correction and treat the difference in each 
as a lcr error. Similarly, for the perturbing cluster 
galaxies, we assumed a standard elliptical spec- 
trum to convert to a B band luminosity. For the 
cluster redshift of z = 0.33, the colour correction 
is 2.06 mag and the k-correction -0.63 mag. 

3.3.1. Modelling procedure 

Since the nearby foreground cluster shears the 
light from the lens galaxy as well as the ring im- 
age, the morphology of the baryons cannot be 



simply fixed as the Sersic profile fit to the ob- 
served lens galaxy surface brightness (see Section 
[2]). In theory we could allow the baryonic mor- 
phology to vary in the minimisation. However, in 
practice, the morphology proves very poorly con- 
strained since the baryons have only a small elon- 
gation (see Section 14. 1|) and only contribute ap- 
proximately one third of the total projected mass 
inside the Einstein radius (including the cluster 
convergence) . 

With this in mind, we adopt the following pro- 
cedure when fitting the lens with the dual compo- 
nent model and pure M/L model: 

1. First fit the gNFW, power-law and SIE mod- 
els to determine the best fit cluster shear. 

2. Using the average of this best fit shear, com- 
pute the elongation and orientation of the 
Sersic profile that would have been observed 
in the absence of the cluster. 

3. Fit the dual component model, fixing the 
baryon component with the Sersic profile pa- 
rameterised in Table [TJ except with the elon- 
gation and orientation computed in step 2. 

We compute an error on the elongation and ori- 
entation estimated in step 3 by combining the er- 
rors from the original Sersic fit, as given in Tabic 
[TJ and the error on the shear determined in step 1. 
This resulting uncertainty is incorporated into our 
modelling using a Monte Carlo simulation which 
randomly samples the orientation and elongation 
for 100 minimisations. The final errors we quote 
combine the spread in minimised model parame- 
ters with the formal errors from the fit and addi- 
tionally include the colour and k-correction error 
for the M/L. 

Having fit each of the five lens models, we re- 
peat the analysis with the best fit model and a 
second source plane. In this case an additional 
parameter, z-i, is introduced for the second source 
plane redshift. 

4. Results 

4.1. Lens reconstruction 

Table [2] compares the \ 2 f° r the five models. 
The best fit to the observed ring is given by the 
dual component model. The pure M/L model is 
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very strongly ruled out. The pure halo, power-law 
and SIE models give a satisfactory \ 2 but perform 
significantly worse than the dual component model 
with Ax 2 = 17.4 for six fewer degrees of freedom 
for the gNFW, A% 2 = 16.2 for 12 fewer degrees of 
freedom for the power-law and A\ 2 — 22.8 for the 
SIE. 

We fitted the gNFW, power-law and SIE mod- 
els first to establish the cluster shear (see Sec- 
tion I4.1.ip . The normalisation of the SIE model 
corresponds to a central velocity dispersion of 
210 ± 10 km s _1 , in agreement with the mea- 
surement of ~ 230 ± 30 km s _1 from the lens 



galaxy spectrum (jSmail et al. 1120071 ). The slope 



of the gNFW model is 1.95 ±0.03, consistent with 
the slope of 2.01+ ° ° 2 averaged over 15 lenses by 



Koopmans et al.l ( 2006^ 1 . Similarly, the best fit 
slope of the power-law model is 2.09 ± 0.04. 

Using the average shear obtained from the 
gNFW and SIE models, we distorted the Sersic 
profile fitted to the observed galaxy light to give 
a resulting elongation of 1.14 ± 0.08 and orienta- 
tion 23.2° ± 4.2°. This sheared Sersic profile was 
then used to fix the morphology of the baryonic 
component for the pure M/L model and the dual 
component model. 

The best fit parameters of the dual compo- 
nent model are given in Table [H The baryonic 
component accounts for 46 ± 11% of the total 
lens mass inside the Einstein radius of 0.77" (5.6 
kpc). This excludes the convergence from the clus- 
ter and its nearby perturbing galaxies which con- 
tribute a total of ~ 20% of the projected mass in- 
side the Einstein radius. The baryonic fraction is 
lower than the fraction 64 ±4% (18% rms scatter) 



Model 


Xmin 


NDoF 


N src 


N par 


dual component 


1912.6 


1896 


278 


10 


SIE 


1935.4 


1888 


289 


7 


pure gNFW 


1930.0 


1890 


285 


9 


power- law 


1928.8 


1884 


291 


9 


pure M/L 


2188.4 


1920 


260 


4 



Table 2: Comparison of the lens models in terms 
of \ 2 a t the minimum. The number of degrees 
of freedom (NDoF) for each model is given by the 
number of image pixels in the elliptical annulus (= 
2184 - see Figure [3]) minus the sum of the number 
of source pixels in the adaptive grid (N src ) and the 
number of parameters in the lens model (N par ). 



found by iGavazzi et al. ( 2007 ) in averaging over 
22 strong lenses, although is within their spread of 
40%— 100%. The halo alignment is consistent with 
the baryons, both positionally and in orientation. 
The elongation of both components is also consis- 
tent, although the errors are relatively large, pre- 
dominantly due to a degeneracy between the halo 
elongation and the shear. We discuss the cluster 
shear further in Section 14.1.1 1 

In Figure [21 we plot the confidence levels 
on the fitted inner halo slope and baryonic 
M/L, having marginalised over all other pa- 
rameters. The best fit slope is a — 1.42+Q 22- 
This is consistent with t he value measured by 
iTreu fc Koopmans I (|2004l ) of l-3±g;2 with their 
baryons + halo model, averaged over three lenses 
although slightly higher than the value measured 
by DW05, a = 0.87±gj$ (their 68% CL) for the 
Einstein ring 0047-2808. The best fit M/L is 
Mq/Lbq which comp ares to the local 



-i 7-1+0.28 
i - ' 1 -0.38 



value for ellipticals measured by iGerhardt et al. 
(|200ll) of 7.3 ± 2.1 M Q /L S0 . To evolve this local 



measurement back to the redshift of o ur len s, we 
use the findings of iKoopmans et al. that 



dlog(M/L B )/dz = -0.69 ± 0.08. This gives a re- 
sult of M/Lb = 2.3 ±0.7 M Q /L B0) slightly higher 
than our measurement but consistent within the 
errors. 

A source of error not included in Figure [2] is the 
uncertainty on the halo scale radius ?'o. However, 
the propagation of this uncertainty is small. We 
find that a 10% change in tq gives rise to a ~ 
1% change in the minimised M/L and a negligible 
change in the minimised slope. This is exactly the 
same dependency found by DW05. 

Finally, turning to the perturbing cluster galax- 



Parameter 



Minimised Value 



Ps 

a 
4' 

eh 

(Ax h ,Ay h ) 


(1.42 ± 0.24) x 10' Mokpc"- 1 
1 42+ - 24 
1.7ll°;l M /L S0 

i i 9 +0.09 
— 0.07 

(-0.03" ± 0.03", 0.06" ± 0.03") 


6 h 


20.3° ± 3.5° 


7 


0.181 ±0.021 


0~.. 


29.6° ±4.1° 


*c 3 


3.6 ± 1.9 M Q /L BQ 



Table 3: Best fit parameters for the dual compo- 
nent model (see Section r3.3j) . la errors are quoted. 



S 



ies, we find a best fit B band M/L of ^> cg = 
3.6 ±1.9 Mq/Lbq- This is in keeping with lo- 
cal ellipticals as well as ellipticals found in z ^0.3 
clusters (e.g., iNataraian fc Springell |2004 ). We 
note that if these cluster members are omitted 
from the modelling then the best fit elongation 
of all models (apart from the pure M/L) has to be 
significantly higher, ~ 1.4 — 1.5 in order to fit 
the observed ring image. In this case, the fit is not 
significantly degraded, but contradicts the concor- 
dant view that isodensity contours have consistent 
or smaller elongatio ns than isophotal contou rs for 



a given system (e.g., Koopmans et al. 2006). 



Cluster shear 

The best fit cluster shear predicted by the 
gNFW, power-law and SIE models was found to be 
respectively (7 = 0.182± 0.017, 7 = 27.6°±3.3°), 
(7 = 0.180 ± 0.018, 7 = 28.1° ± 3.4°) and 
(7 = 0.176 ± 0.014, 7 = 29.1° ± 2.9°). 

As well as being self-consistent, this shear is in 
agreement with what would be expected from an 
isothermal cluster mass profile: Firstly, the tan- 
gential shear angle 7 is consistent with the di- 
rection of the cluster centre at 24°. Secondly, the 
magnitude is what would be expected given the 




Fig. 2. — One-parameter confidence regions (1, 2, 
3 & 4<r) on the halo inner slope, 7, and M/Ls© of 
the baryons, for the dua l component model (see 



Treu fc Koopmans I 12004 iDve fc WarrerTI 120051 



for comparison). 



cluster proximity as the following approximate cal- 
culation demonstrates: The ring lies ~ 1.8 times 
further away from the cluster centre than the out- 
ermost (and faintest) arc. Making the reasonable 
assumption that this arc is an image of a z > 2 
source and since the deflection angle varies very 
slowly with rcdshift beyond z ~ 2, this gives a 
good approximation to the critical radius for a 
source at infinity. Treating the cluster mass profile 
as an isothermal sphere gives a predicted shear of 
7 = 0.28 for a z — > 00 source at r = 1.8r cr it- How- 
ever, a source at the lens redshift of z = 0.73 expe- 
riences a shear of 0.63 times less (for our assumed 
cosmology) resulting in a shear of 7 = 0.175. 

The average of the shear from the gNFW, 
power-law SIE models was used to distort the Ser- 
sic profile fitted to the observed galaxy light. With 
this distorted profile, the dual component model 
gives a best fit cluster shear of 7 = 0.181 ± 0.021, 
7 = 29.6° ± 4.1°. The error budget here in- 
cludes the uncertainty on the distorted Sersic pro- 
file. This shear agrees very well with the shear ob- 
tained from the gNFW, power-law and SIE mod- 
els, an indication that the lack of baryons in these 
models did not bias the shear they predict. 

The ability to predict the large scale cluster 
shear from the small area of sky covered by the 
ring is quite remarkable. Although the shear is 
fairly degenerate with the lens elongation, this 
degeneracy is significantly higher in strong lens 
systems with only a point-like source, to the ex- 
tent that measurement of external shear in these 
systems is impossible. This is a testament to 
the stronger constraints provided by an extended 
source. 

4.2. Source reconstruction 

Figure [3] shows the reconstructed source ob- 
tained from the best fit dual component model. 
The top row corresponds to the unregularised so- 
lution and the middle row the regularised case. In 
the bottom row we plot the observed image along- 
side the image of the unregularised source and the 
significance of the residuals. The residuals are 
within ±2.5(7. Considering some of the pixels in 
the observed image reach a significance of ~ 40er, 
this demonstrates the quality of the fit. Never- 
theless, we have investigated the cause of these 
residuals, i.e., whether they are due to an overly- 
coarse source pixelisation, the parameterisation of 
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Fig. 3. — The reconstructed source from the best fit dual component model. Reading from left to right along 
each row; Top row - the unrcgularised reconstructed source, the standard errors map for the reconstructed 
source and the significance map; Middle row - same as top row but for the regularised source; Bottom row - 
the masked observed ring, the image of the unrcgularised source and the significance of the residuals. As can 
be seen, the lens model gives a very good description of the observed ring morphology. The reconstructed 
source shows two peaks. In the top left panel, the dashed squares delineate the two source planes used in 
the dual-source reconstruction (labeled 1,2). The source plane caustic is indicated in the top two rows. In 
the source plane, 0.1" corresponds to 800pc. 
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the lens model or the PSF used to smear images 
of lensed source pixels. The source pixelisation is 
ruled out by the fact that the residuals are barely 
changed when we use the image of the regularised 
source. We also repeated the minimisation with 
different PSFs, including a PSF extracted from 
stars in the field, but the best fit was obtained 
with the TinyTim model PSF (see Section [2]). We 
conclude therefore that the residuals would prob- 
ably be lessened by a small modification in the 
parameterisation of the lens model. This would 
require an exhaustive search through different pa- 
rameterisations with only a small return and no 
guarantee of finding a unique solution, hence we 
leave this for possible future work. 

The total magnification of the system (i.e., to- 
tal ring flux divided by total source flux) is ~ 25. 
There are clearly two peaks in the source surface 
brightness distribution. To demonstrate the con- 
tribution each peak makes to the observed Ein- 
stein ring image, Figure [4] shows the image of each 
peak lensed separately with the best fit dual com- 
ponent model. The dominant ring structure is due 
to the brighter peak mostly contained within the 
caustic. The fainter peak lying nearer the top of 




D 0.2 0.4 0.6 0.2 0.4 0.6 

X (orcsec) X (arcsec) 




1 2 3 1 2 3 



X (arcsec) X (arcsec) 

Fig. 4. — The contribution of each reconstructed 
source to the observed ring image. The top row 
shows the two sources and the bottom row their 
respective images. 



the plane and outside of the caustic is imaged into 
the westerly and easterly extensions. Although 
not plotted, the double-peaked nature of the re- 
constructed source obtained with the other models 
is reproduced in every case. 

4-- 2.1. Dual source plane analysis 

We repeated the minimisation using the dual 
component model but with two separate source 
planes as indicated in the top left panel of Figure 
El Source plane 1 containing the brighter source 
peak was held fixed at the redshift z = 3.07 but 
the redshift of source plane 2, Z2, was allowed to 
vary in the minimisation. 

Our results give a value of Z2 = 3.13^Q4g with 
an overall fit that is not an improvement on the 
single source plane case. The remaining minimised 
parameters of the model differ only negligibly to 
those in Table for the single source plane case. 
Although the errors are large, the redshift of the 
second source peak is consistent with that of the 
brighter primary peak, suggesting that these are 
two features within the same galaxy. 

5. Summary and Discussion 

Out of the five lens models assessed in this 
work, the Einstein ring LBG J2135-0102 is best 
fit with our dual component model comprising a 
gNFW dark matter halo that hosts a baryonic 
component following the lens galaxy light. Of 
the remaining four single component models, the 
gNFW, power-law and SIE models give an accept- 
able fit but the pure M/L model is very strongly 
ruled out. Between the gNFW, power-law and 
SIE, the gNFW gives a better fit with Ax 2 = 1.2 
for six more degrees of freedom compared to the 
power-law and A\ 2 = —5.4 for two more degrees 
of freedom compared to the SIE. However, all 
three models give a significantly worse fit than the 
dual component model which has A% 2 = —17.4 
for six more degrees of freedom compared to the 
gNFW model. 

The dual component model predicts a projected 
baryonic contribution of (46 ± 11)% interior to the 
Einstein radius of 0.77". This is at the low end 
(but consistent with) the spread i n baryonic frac- 
tion o f 40% - 100% measured by iGavazzi et al. 
( 20071 ) over 22 strong lens systems. We measure 
a M/L ratio of Ull^g M Q /L BQ for the bary- 
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onic component, in agreement with the value for 
local ellipticals once th e evolution measured by 



Koopmans et al. 



( 20061 ) is taken into considera- 
tion. 

We have modelled the effect of the z — 0.33 
foreground cluster MACS J2135. 2-0102 on our lens 
solution, both in terms of the large scale shear 
and convergence and also perturbations to the lo- 
cal potential from nearby cluster members. The 
modelling accurately predicts the cluster shear. In 
strong lens systems with only a point-like source, 
this is impossible due to a strong degeneracy be- 
tween the shear and lens elongation. Although a 
degeneracy exists between the cluster shear and 
elongation in our modelling, the extra constraints 
provided by the ring image greatly reduce it. This 
demonstrates the advantage extended source sys- 
tems give over point source systems. 

In addition to shearing the ring image, the clus- 
ter shears the lens galaxy light. Taking this shear 
into account, we find that the halo and bary- 
onic component are well aligned (i.e., in terms of 
the centroid and orientation of the major axis) 
within uncertainties. Failure to incorporate the 
local perturbing cluster members does not signif- 
icantly degrade the fit but results in a substan- 
tially higher halo elongation (e/, ~ 1.4 — 1.5) than 
that of the intrinsic lens galaxy light. This per- 
turbation is therefore required to prevent con- 
tradiction with the accepted view that a given 
galaxy's isophotes should be consistent with or 
more elongated than it s isodensity contours (e.g., 



Koopmans et al. 2006h 



The best fit inner slope of the total mass is 
a = 1.95 ± 0.03, given by the gNFW model 
or a — 2.09 ± 0.04 given by the power-law 
model. This is in keeping with the findings of 
seve ral strong lens st u dies to date, for exam- 



ple Koopmans et al. ( 20061 ) who measured a 
slope of 2.0ltnn^ averaged over 15 lenses and 



Rusin. Kochanek fe Keeton I (|2003f ) who measured 
a slope of 2.07 ±0.13 averaged over 22 lenses. The 
result lends further evidence towards the proposi- 
tion that the slope has little or no evolution out 



to z ~ 1 ( Koopmans et al. 2006f ). 



With the dual component model, the best fit in- 
ner slope of the halo component is a = 1A2~Tq22- 
This is consistent with the value of 1.3 1°'^ av- 



(|2004l ), although slightly higher than the value 
a = 0.87l°;27 measured for the z = 0.485 lens 
0047-2808 by DW05. To address whether this is 
consistent with pure CDM simulations, the ef- 
fect of halo contraction by t he condensation of 
baryo ns must be considered 



Blumenthal et al. 



( 19861 ) originally suggested the model of ad i abatic 
contraction. More recently. iGnedin et al. ( 20041 ) 
showed using high resolution CDM + baryon sim- 
ulations of halos that the adiabatic contraction 
model over-predicts the in crease in central dark 



matt er density (see also ISelwood fc McGaugh 



2005). They provide analytical fitting functions 



eraged over three lenses by iTreu fc Koopmans 



to describe how a NFW profile is contracted by a 
collapsed baryonic profile of arbitrary inner slope. 
The effective inner slope of the volume mass den- 
sity profile corresponding to our baryonic Sersic 
law is ~ 1.7 which, according to these fitting func- 
tions, would contract a NFW profile into a profile 
with a slope of ~ 1.5. The slope determined in 
our analysis for LBG J2 135-0 102 therefore corre- 
sponds to an uncontracted slope of ~ 1, consistent 
with current CDM simulations. 

This paper presents detailed modelling of the 
inner halo mass profile of only one Einstein ring 
system. To date, there are a further ~ 25 such 
systems imaged with the HST for which detailed 
dual component models have not yet been fully 
published. Analysis of these remaining systems is 
crucial to provide a better measure of the mean 
and dispersion of the inner halo slope thus greatly 
improving current strong lensing constraints on 
the CDM model. 

Finally, considering the reconstructed source, 
we can conclude that the two peaks seen in the sur- 
face brightness distribution are most likely struc- 
ture within the same galaxy. The magnification of 
the ring allows the source to be reconstructed with 
a pixel scale of ~ 1/2 that of the observed image at 
the same signal-to-noise and with approximately 
non-covariant pixels. This gain in spatial resolu- 
tion means that integral field spectroscopy in the 
optical and more importantly with laser adaptive 
optics in the near infra-red offers the opportunity 
to spatially resolve the star-formation, kinematic 
and chemical properties of individual Hn regions 
within t his galaxy on scales of ~ 200p c (for exam- 
ple, see ISwinbank et al. I [20061 120071) . This level 
of science is already on a par with that which will 
be delivered by the next generation of Extremely 
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Large Telescopes. 
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